Decoherence, entanglement decay, and equilibration produced by chaotic 

environments. 
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We investigate decoherence in quantum systems coupled via dephasing-type interactions to an ar- 
bitrary environment with chaotic underlying classical dynamics. The coherences of the reduced state 
of the central system written in the preferential energy eigenbasis are quantum Loschmidt echoes, 
which in the strong coupling regime are characterized at long times scales by fluctuations around 
a constant mean value. We show that due to the chaotic dynamics of the environment, the mean 
value and the width of the Loschmidt echo fluctuations are inversely proportional to the quantity 
we define as the effective Hilbert space dimension of the environment, which in general is smaller 
than the dimension of the entire available Hilbert space. Nevertheless, in the semiclassical regime 
this effective Hilbert space dimension is in general large, in which case even a chaotic environment 
with few degrees of freedom produces decoherence without revivals. Moreover we show that in this 
regime the environment leads the central system to equilibrate to the time average of its reduced 
density matrix, which corresponds to a diagonal state in the preferential energy eigenbasis. For the 
case of two uncoupled, initially entangled central systems that interact with identical local quantum 
environments with chaotic underlying classical dynamics, we show that in the semiclassical limit the 
equilibration state is arbitrarily close to a separable state. We confirm our results with numerical 
simulations in which the environment is modeled by the quantum kicked rotor in the chaotic regime. 

PACS numbers: 



I. INTRODUCTION 

In open quantum systems the interaction between a 
system and its environment may result in the well known 
phenomenon of decoherence 0, While the central 
system becomes ever more entangled with the environ- 
ment, quantum information initially present in the re- 
duced state of the system may be lost to the environment. 
The most transparent example of how the irreversible 
loss of quantum information to the environment leads 
to decoherence is in the case of dephasing-type system- 
environment interactions (or so-called measurement-type 
interactions^), i.e., when the interaction Hamiltonian 
commutes with the free system Hamiltonian. In this case, 
by tracing over environmental degrees of freedom, one 
may observe the irreversible decay of the quantum co- 
herences of the central system's reduced density matrix 
written in the preferential basis of the free Hamiltonian 
eigenstates, while the populations are conserved. 

The traditional approach to this problem points to the 
need for an infinite number of environmental degrees of 
freedom in order for the decoherence process to occur Q ■ 
The Caldeira-Leggett model [|| is the most renowned en- 
vironment model of this type. However, recently there 
has been an increasing interest in understanding if and 
how environments with few degrees of freedom can pro- 
duce decoherence in a quantum system of interest 0- 
[l8j and if decoherence can be produced by few internal 
degrees of freedom [ljj [20[. One motivation for this is 
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the relevance of decoherence to quantum computation 
and quantum information tasks, where often the interac- 
tion of the system of interest with a global environment 
(usually composed of many degrees of freedom) is well 
screened and therefore the interaction with a "near" en- 
vironment composed of few degrees of freedom involved 
in the control of quantum operations may be the most 
relevant [U, [22| . The other motivation is the study of 
the emergence of classicality in quantum systems, where 
decoherence plays a central role. 

The study of environments with few degrees of freedom 
is in general related to the role of the chaotic dynamics 
of their classical counterpart [H|[l6j]. The con- 

nection between environment-induced decoherence and 
the so called fidelity decay (or quantum Loschmidt-echo 
decay J47|) brought great insight into the study of de- 
coherence by chaotic environments [l|, [Tfjl mum. 
The sensitivity of chaotic environments to the pertur- 
bations produced by the interaction with a central sys- 
tem was conjectured in 23] to be related to their abil- 
ity to rapidly produce decoherence. A precise connec- 
tion between quantum Loschmidt echo and decoherence 
was first made in [25| in the case of dephasing-type 
system-environment interactions, where the off-diagonal 
elements of the system's reduced matrix written in the 
preferential basis are different Loschmidt echoes in the 
environmental degrees of freedom. This connection does 
not neglect the free evolution of the central system. A 
similar approach was developed in [l3j , but for a generic 
small coupling to the environment. In that article the 
authors show, using perturbative approximations, that 
the off-diagonal elements of the reduced matrix of the 
system in its energy eigenstates basis are proportional to 
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fidelities. 

It is common to associate the decoherence time with 
the characteristic short-time scale for the decay of the 
fidelity [IH, [23|. However, for long times the fidelity in 
general can suffer fluctuations that can be large and may 
lead to important revivals of the coherences. Indeed, if 
the quantum environment has chaotic underlying classi- 
cal dynamics, at long-time scales the coherences of the 
central system's reduced density matrix, written in its 
preferential eigenbasis, fluctuate around a constant mean 
value. For a generic central system strongly coupled by 
dephasing-type interactions to an arbitrary chaotic en- 
vironment (which may have few degrees of freedom) , we 
show that the time-average and the width of these fluctu- 
ations are inversely proportional to the effective dimen- 
sion of the environment's Hilbert space, which can always 
be defined (even in the case in which the total available 
Hilbert space of the environment has infinite dimension) 
. Thus, for chaotic environments, the decoherence oc- 
curs in the semiclassical limit, i.e., h c ff = h/S — > (S is 
a typical action of the environment), where typically the 
effective Hilbert space dimension of the environment is 
large. 

Our results have a direct application to the problem 
of equilibration, a key process in the understanding of 
thermal equilibration in quantum systems (30l - [3^ |. A 
central system equilibrates if its initial state evolves to- 
ward some particular state, in general mixed, and re- 
mains in that state, or close to it, for all times. Here 
we show that if the environment Hamiltonian has un- 
derlying classical chaotic dynamics, and for a dephasing 
system-environment interaction, the central system equi- 
librates to a totally decohered mixed state in the semi- 
classical regime, independent of the initial state of the 
system and for generic initial states of the environment 
(initially decoupled from the system of interest). 

The study of entanglement decay due to the action of 
environments is a central problem in quantum compu- 
tation and quantum-information processing. In this re- 
gards it is important to know how chaotic environments 
with few degrees of freedom can produce ent ang lement 
decay in the system of interest d, [111, [HI, EE EI[. Here 
we show that an initially entangled state of two non- 
interacting central systems coupled to equivalent local 
quantum environments with chaotic classical dynamics 
equilibrates to a state that is arbitrarily close to a sepa- 
rate state in the semiclassical limit. 

The paper is organized as follows: In Sec|H]we present 
the echo dynamics approach to decoherence. We then ob- 
tain the long-time behavior of decoherence functions for 
chaotic environments in Sec lIIII and go on to discuss de- 
coherence and equilibration of a one party central system 
in Sec lIVI Disentanglement and equilibration of bipartite 
systems produced by chaotic environments are discussed 
in SecfVland in Sec lVII we describe our numerical simula- 
tions and analyze their results. We finish with concluding 
remarks in Sec lVIII 



II. ECHO DYNAMICS APPROACH TO 
DECOHERENCE 

We consider an arbitrary system interacting with a 
generic dephasing environment: 



H = H c <8> l e + gS ® V + l c ® H e 



(1) 



where H c and H e are, respectively, the Hamiltonians of 
the central system and the environment. The dephasing 
interaction is Hi = gS ®V, where S acts on the system 
degrees of freedom and [H C ,S] = 0. The operator V 
acts on environmental degrees of freedom, and g is the 
coupling strength. The global initial state is a product 
state: 

p ce (0)=p(0)®o>(0)=^A„ TO |n)(m|®(2i(0) , (2) 



where w(0) is the initial reduced state of the environment, 
and we have expanded the initial state of the central sys- 
tem p(0) in the preferential basis of common eigenstates 
of H c and S, which we assume to have discrete and non- 
degenerate spectra: 



H c \n) = e„|n), 
S\n) = s n \n). 



(3) 



The evolved reduced density matrix of the central sys- 
tem is obtained by tracing over the environmental degrees 
of freedom 

p(t) = Y J A nm e-^-^ t ' h Tr e [ulU n w(0)\ \n)(m\, 

n,m 

where t/f+ e , rl ( m ) = C^ n (m) are conditional effective evolu- 
tion operators associated with the Hamiltonians 



TT JT 



e H~ .9 s ra(m) 



V, 



(5) 



which act exclusively on environmental degrees of free- 
dom. Hence, the central system's diagonal matrix ele- 
ments in the preferred basis {\n)} arc constant, while 
the time evolution of off-diagonal matrix elements is con- 
trolled by the decoherence functions @ 



^nm (^) 



Tr e U m {-t)U n (t)u(0) ,n^m. (6) 



Following essentially the approach in [25[ (see also [lj, 
El, H3 ) , we can associate the decay of these decoherence 
functions with the decay of quantum Loschmidt echoes 
in the Hilbert space of the environment. Indeed, the 
amplitudes of the decoherence functions © are the so- 
called allegiance amplitudes [19 ]. 



fnm{t) = Tr e M nm (t)Cj{0) 



(7) 



where we have introduced the echo operator that acts in 
the Hilbert space of the environment, 



M nm {t) = U m (-t)U n (t). 



(8) 
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In this echo dynamics, H m (J5)) plays the role of the un- 
perturbed Hamiltonian, and H n ([5]), rewritten as H n = 
H m + tnmV i is the perturbed Hamiltonian with pertur- 
bation amplitude given by e nm = g(s n — s m ). 

For pure initial states w(0) the allegiance amplitudes 
in Eq.([7|) reduce to the so-called fidelity amplitudes in- 
troduced by Peres [H]. It is important to note that the 
decay of an off-diagonal element in a given column of the 
system's reduced density matrix is determined by an echo 
operator composed of an unperturbed evolution operator 
U m , which is common to all the elements in that column 
and a perturbed operator U n . The perturbation strength 
e nm increases as we move away from the diagonal. Thus, 
all echo dynamics associated to a given column m are in 
the strong perturbation regime if and only if e m +i. m V 
represents a strong perturbation. 

In general, the time evolution of the F nm (t) = 
\fnm(t)\ 2 depends on the specific system-environment 
coupling, on the environment Hamiltonian H e , and on 
the properties of the initial state w(0). Complete deco- 
herence is said to have occurred if after a typical time 
scale , called the decoherence time, the system's reduced 
density matrix becomes diagonal in the preferred basis. 
Nevertheless, we can also say that the central system suf- 
fers decoherence or partial decoherence if in the long-time 
regime, F nm (t) decreases to very small values and re- 
mains small for times much longer than the decoherence 
time (the central system never regains its coherences). 
This means that in order to produce true decoherence 
on another system, the environment echo dynamics must 
decrease significantly in the long-time regime and must 
not present revivals thereafter. 

A very simple example of a rapid initial decrease of 
F n m(t) which presents revivals in the echo dynamics con- 
sists of an arbitrary central system coupled in the form 
Hi = H c ® Y^jLiidjb] + 9jbj) to a non-chaotic environ- 
ment composed of N bosonic modes described by a set of 
harmonic oscillators, H e = Y^jLi ttojbjbj. When the ini- 
tial state of the environment is w(0) = |0) (0| ® . . .<8 |0) (0| 
(where |0) is the vacuum state), and for a constant spec- 
tral density of environment modes with a finite cut-off, 
it is straightforward to show that the system coherence 
amplitudes will oscillate in time. So, at finite times the 
central system regains its lost coherence. In order for de- 
coherence to occur, without revivals, one must perform 
the continuous limit of environment modes. In this case 
the decoherence functions are given by exponentially de- 
creasing functions with no revivals. Therefore, in this 
simple example, the environment can produce true deco- 
herence on the condition that it has a very large number 
of degrees of freedom. 

In the case of an environment with few degrees of free- 
dom, we will show that a very similar situation arises 
when the echo dynamics in the environment is associated 
with a Hamiltonian which has chaotic underlying classi- 
cal dynamics. In this case decoherence (with no revivals) 
is produced when the effective dimension of the environ- 



ment's Hilbert space is very large, a condition that gener- 
ically can be satisfied in the semi-classical regime (small 
effective Planck constant), independent of the system- 
environment initial state. 



III. LONG-TIME BEHAVIOR OF 
DECOHERENCE FUNCTIONS F nm (t) FOR 
CHAOTIC ENVIRONMENTS 

In this section we determine the long-time behav- 
ior of the decoherence functions F nm {t) in the case of 
chaotic environments. When the classical dynamics as- 
sociated with the free environment Hamiltonian H e is 
fully chaotic, one can generally expect that the classical 
dynamics associated with the effective Hamiltonian H m 
[Eq.©] is also fully chaotic. In this case, the behavior 
in time of the decoherence functions F nm (t), determined 
by the echo operator in Eq.©, consists essentially in an 
exponential decay with different decay rates depending 
on the perturbation regime [H, H3- However, when the 
Hilbert space of the environment is finite, the decoher- 
ence functions will not decay to zero even at long-time 
scales. Indeed, the discrete spectrum of the evolution 
operators in the echo operator in Eq.© causes F nm to 
fluctuate around the time-average value (F nm ) (27l. [33|. 
where 

(.) = hm - f mdt. (9) 

Nevertheless, even if the Hilbert space of the environ- 
ment is infinite, the effective Hilbert space covered by 
the evolved state of the environment is always limited, 
in which case (F nm ) also has a finite value. A straight- 
forward example is when the average available energy is 
finite. In this case, the state of the environment is con- 
strained to the subspace composed of energy eigenstates 
with less energy than this available amount. 

With the assumption that the quantum echo dynam- 
ics in the environmental degrees of freedom is associated 
with fully chaotic underlying classical dynamics, we ex- 
tend the arguments in [27| to obtain a relation between 
the time average of the decoherence functions and the 
effective dimension of the environment's Hilbert space 
in the strong perturbation regime. Indeed, by expand- 
ing the initial state of the environment in the eigenbasis 
{|?/>z)} of the unperturbed evolution operator or the 
eigenbasis U 4>) I of the perturbed evolution operator E/t 
one obtains [481 ] : 

N m N m 

*(°) = £X>*( )hfoW*l= ( 10 ) 

1=1 fc=l 

= E E c '*(°W'>$*i' ( n ) 
i=i k=i 

where N m and N nm represent the number of eigenstates 
of the evolution operators which significantly contribute 
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to the sums, i.e., uiik{0) « for I, k > N m and uiikifi) ~ 
for l,k > N nm . At any future time, the reduced den- 
sity matrix of the environment ui(t) can be spanned by 
iVefi linearly independent orthogonal eigenvectors {IV'z)} 
or {IV 1 /)}, where 



N. 



e ff = max(iV m , AT nm ). 

nm 



(12) 



Therefore, N e g can be interpreted as the effective dimen- 
sion of the environment's Hilbert space for all times. In 
Appendix [A] we show that due to the classically chaotic 
underlying dynamics of the environment, the strong cou- 
pling regime to the central system (and for N e s 3> 1), 
it is possible to relate the time average of the system 
dccohcrcncc functions to: 



(Fnm) = ^ = C 



ff ' 



(13) 



where C is a constant that depends on a)(0) (see Ap- 
pendix [A} and C is related to C once we know the specific 
relationship between the effective Planck constant S e g 
and iVefi that can involve an exponent 7. For example, 
in the case of environments with an autonomous Hamil- 
tonian H e and time-independent couplings to the central 
system, the dynamics associated with the unperturbed 
Hamiltonians H m conserve energy. In this case, follow- 
ing the well known semiclassical prescription (i.e., Weyl 's 
rule [Hj]) one obtains N c g ~ N m fts i/E m /(27rft ff) 7 , where 
i>E m is the volume inside the phase-space surface of con- 
stant energy E m ~ Ti[H m u}(0)], and 7 is the number of 
degrees of freedom of the environment. 

The description of the long-time behavior of the deco- 
herence functions F nm (t) in the case of chaotic environ- 
ments is completed with the analysis of the width of the 
fluctuation around the mean value. In Appendix [A] we 
obtain, 



AF nm = y/{F* m ) - (F, 
G 



N, 



eff 



(14) 



where G is a constant that depends on o>(0) (see Ap- 
pendix [AJ and G is related to G once we know the spe- 
cific relationship between h e a and N e g. We point out 
that to arrive at relations (TIB")) and (fT4")) . the approxi- 
mation N c ff 1 must be made. With these conditions 
fulfilled we can say that the above results only depend 
on the initial state of the environment ui(0) through the 
values of C and G. 



IV. DECOHERENCE AND EQUILIBRATION 
OF A CENTRAL SYSTEM BY A CHAOTIC 
ENVIRONMENT 

Here we show that the results of Sec. IIIII imply that 
in the semiclassical regime (H e s < 1) a central system 
with a dephasing-type coupling to an environment with 



fully chaotic dynamics equilibrates to a state given by 
the time- averaged state of the central system, i.e., 



Pe- 



rn) 



n)(n\, 



(15) 



where pit) is given in Eq.@ [see Appendix iBlfor a deriva- 
tion of Eq. (fl"5"j) ]. This occurs independently of the initial 
state of the environment and of the central system pro- 
vided that the initial state of the composite system is 
uncorrelated (as in Eq.([2])). 

Generically, there is a back and forth of information 
flow between the central system and the environment as 
the reduced state of the system pit) fluctuates around 
the totally decohered state (pit)) in Eq. (|T5|) . As a 
consequence of this non-Markovian character of the re- 
duced dynamics of the central system, the degree of dis- 
tinguishability between the evolved state pit) and (pit)) 
oscillates. In order to quantify this fluctuations we use 
the trace distance 341: 



D(p(t),p eq ) = iTr cV /[/5(t)-p cq ] 



(16) 



which is a measure of the distinguishability between the 
two quantum states [3(|. An upper bound to the trace 
distance is given by the Hilbert-Schmidt distance Dhs 
51, 



DW),P, 



<■-</ 



< 



fN, 



c -D HS (p(t),p eq ) (17) 



IJ2 l^nm| 2F nmO) (18) 



where D HS ip(t), p cq ) = y/Tr c [(p(t) - p eq ) 2 } and N c is 
the effective dimension of the Hilbert space of the cen- 
tral system, which is given by the number of essentially 
non-zero terms in the expansion of the central system's 
initial state p(0) in the eigenbasis of H c . We calculate 
the Hilbert-Schmidt distance in Eq. (fT5)) from the expres- 
sions for p(t) (01, (p(t)) (rr5"|) . Using the concavity of the 
square-root function we can estimate an upper bound for 
the time-averaged fluctuations: 



(D( P (t),p eq )) <h, 



off o 




(19) 



where we use the result in Eq. (fT"3"]) for (F nm (t)). Ac- 
cording to this result, the fluctuations, measured by 
(D(p(t), Peq)), tend to zero in the semiclassical limit 
(H e g — > 0). This means that in this limit the evolved 
state pit) becomes indistinguishable from the totally un- 
correlated state in Eq. (fT"5"| . In the semiclassical regime, 
the entanglement between the central system and the en- 
vironment reaches its maximum possible value (which is 
determined by /5(0)). 

We can see this by noting that if the initial total com- 
posite state is /3 ce (0) = p(0) <8>cD(0), and for a dephasing- 
type coupling of the central system to the environment, 
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the purity of the reduced evolved density matrix has the 
lower bound, 

Tr c [p(i) 2 ] = Yl K™| 2 + E \ A nn\ 2 F nm (t) 

n ri^ra 

> n inv (p(o)) > (20) 

where lZ- mv (p(0)) = Tr[diag/5(0)] = |A„„| 2 is the gen- 
eralized inverse participation ratio [27l ] of the initial re- 
duced density matrix in the energy eigenbasis of the cen- 
tral system. The maximum entanglement between the 
central system and the environment is attained when the 
lower bound in Eq. (|20j) is reached. By inserting the re- 
sult of Eq. (fl"3]) into Eq.(f20|) and taking the time-average, 
one obtains 



where f„ im ,{t) = Tr e 



M„ 



(t)p ei (Q) is the allegiance 



2l\ Seff->0 



^inv(/5(0)) , 



(21) 



which means that maximum entanglement is expected in 
the scmiclassical regime. 



V. DISENTANGLEMENT BY CHAOTIC 
ENVIRONMENTS 

We now consider two noninteracting central systems, 
H Ci (i = 1,2), each coupled to a local environment H ei 
through a dephasing-type coupling, 



H = H Cl + H C2 + H h + H h + H ei + H e2 



(22) 



where Hi i = giSi <E) %, with [H Ci ,Si] = 0. For sim- 
plicity we consider the two local environments to have 
Hamiltonians, H Cl and H £2 , with the same functional 
form. Thus, classically this corresponds to two different 
chaotic systems with the same dynamics. We stress here 
that the occurrence of local environments is the most 
common situation in proposals of quantum computation 
and quantum information processing. In this framework 
the two central systems in Eq. ((22|) can be considered as 
noninteracting qubits in the time interval separating the 
action of two consecutive conditional logical gates that 
involve an interaction between them. Assuming that the 
two system Hamiltonians H Cl and H C2 have discrete and 
non-degenerate spectra, we denote the basis of 

common eigenstates of H Ci and Sf. H Ci \rii} = e ni \rii) 
and Si\rii) = s ni \n.i). The total initial state is the tensor 
product: p ce (0) = u)(0) ® p(0) <g> w(0) where p = p Cl +c 2 is 
the entangled initial state of the bipartite central system 
and d)(0) is the initial state of each local environment. 

After tracing over the environmental degrees of free- 
dom, one obtains the reduced density matrix at time t 
for the central bipartite system: 

p(t) = Tr ei Tr e2 [p ce (t)] 

EA ^(e ni -e mi +e n2 -e m2 )t 
n ?iimiii2m2 ° 

ni ,mi ,n 2 ,rn,2 

x fn imi ( t )fn 2 7n 2 (t)\n 1 )(m 1 \ <g) \n 2 )(m 2 \, (23) 



amplitude introduced in Eq.©, and M nimi are the echo 
operators in the Hilbert space of the local environments 
[see Eq.©]. 

Note that when F niTni = |/„ im J 2 = 0, for all rii ^ nii, 
the initial entangled state p(0) becomes separable. But 
for F nimi ^ the two subsystems are not necessarily 
entangled. In order to illustrate how the allegiance am- 
plitudes fmrni can control the degree of entanglement 
of the central systems, we consider a simple case of two 
qubits with the free system Hamiltonians in Eq. (|22[) given 
by H Ci = huia Zi /2 (i — 1,2), where a Zi are the Pauli 
z operators whose eigenstates are {|0), |1)} ("up" and 
"down" states, respectively). If we start with the entan- 
gled two-qubit state \ip c ) = (|0)|+) + |1)|-)) /2 (where 
|±) = (|0) ± |1)) /\/2), the evolved reduced density ma- 
trix in Eq.dMl) in the basis {|00)|01)|10)|11)} is, 



/ 1 /oi(*)e- 



foi(t)e~ iu 

1/oiWI 2 
1 



V 



h.c. 



-/o 2 i(i)e- 42 ^ 
-/oi(*)e-^ 
-foi(t)e-^ 
1 



(24) 



A simple way to measure the entanglement of this state is 
to calculate its negativity [13], Ne{p) = — Aj)/2, 

where Xj are the eigenvalues of the partial transpose of p. 
For < F 01 (t) < 1, one obtains M E {p) = - (/?)] /8 
where /3 = -F 01 (t) - 2^F 01 (t) + 1 . This means that 
when Foi(t) decreases from 1, the entanglement between 
the two qubits decreases until this function reaches a crit- 
ical value F^ r} = \] -\ + ^2 (when N E {p) = 0) and 
from then on the two qubits are no longer entangled to 
each other. This is a simple example of the so called en- 
tanglement sudden death [38[ by a dephasing reservoir. 

If Foi(t) were a monotonously decreasing function of 
time, we could be sure that from the time in which 
Foi(t) < F^ the two qubits would be disentangled. 
But the disentanglement is also guaranteed even if, af- 
ter some time, -foi(t) starts to fluctuate around a mean 
value (Fqi) < Fq^ with the condition that the width 
AFoi(t) of this fluctuation is sufficiently small. As shown 
in Sec. (|IIl[) . this situation is always satisfied for chaotic 
environments with few degrees of freedom where the 
mean values (F nimi ) and the magnitude of the fluctu- 
ations AF niTni (t) go to zero in the limit of small effective 
Planck constant of the environment. 

We can outline the disentanglement of the two central 
systems due to the action of the local chaotic environ- 
ments if we assume that they are identical, i.e., H Cl and 
H C2 with identical functional form, and with a spectrum 
of eigenenergies with non-degenerate gaps [ijj. In this 
case it is straight forward to calculate the time-average of 
the evolved reduced state of the bipartite central system 
if the initial state has the form /5 ce (0) = w(0)(g)/5(0)<g>a>(0) 
(p = p Cl +c 2 is the entangled initial state of the bipartite 
central system), which unlike the time-averaged reduced 
state in Eq (fT5j) is not necessarily a separable state. Nev- 



G 



crthcless, it can be written as: 

(p(t)) = (p(t))d + 0(h2 s ), (25) 

where 

(p(t))d = EE j4 ™inir l 2« 2 l"-i)( ri il ® \n2){n 2 \ , (26) 

111 "2 

is a disentangled state, and 



1 (i 7, n 1 n 2 (t))|ni)(n 2 | <E> |n 2 )(ni 



(27) 

Note that using Eq. (fT3|) we can write 0(h2 s ) = fi-eff^ 

where Tr[0] = 0. Following the same reason as in 
Section IIVI we can estimate the time-averaged fluctua- 
tions of the evolved state p(t) around (p(t)) and we ob- 
tain: (D(p(t), p cq )) < 0(h2 s ). Thus, we see that in 
the semiclassical regime (p(t)) is the equilibrium state 
of the bipartite central system. It is clear from Eq. (f2"5")) 
that if the equilibrium state is not disentangled it is ar- 
bitrarily close to a separable state in the sense that a 
small perturbation of order 0(h c g ) is enough to separate 
the state. Indeed, if the perturbation is given for ex- 
ample by the completely positive trace preserving map: 

(p) -> (l-<;)(p)+s((p}-6) = {p) d + 6-S, it is enough 
to apply a perturbation of the order <; 
the disentangled state (p)d- 



fQa to obtain 



VI. NUMERICAL SIMULATIONS 

In order to confirm the analytical results in sections 
and IIV| we perform numerical simulations in which 
the environment is modeled by the periodically kicked 
Hamiltonian 



H e = + V Q COs(fco^) J>/ T - n ) 



2M 



(28) 



n=0 



where [X, P] = ih. We consider two forms of dephasing- 
type coupling (see Eq.flT])) to a generic central system. 
The first is a linear coupling via the environment operator 

V = P and the second is via the kicked coupling operator 

V = SV cos(k X)J2^ =0 (t/T — n). Because the results 
we found for the two couplings are equivalent we only 
show those for the linear coupling. 

It is more convenient to work with the dimensionlcss 
coordinate 6 = k$X and momentum p = koTP/M = 
h c ffP/hko that satisfy the commutation relation [9,p] = 
ih c g where the effective Planck constant is h c g = 
kyhT/M. Using these dimensionless variables and per- 
forming the energy transformation (h c gT/h)H — > H to 
the total Hamiltonian H in Eq.((T]), the environment is 
represented by the kicked rotor Hamiltonian (40| . 



KR 



v_ 

2 



(29) 



where the new dimensionless kicking amplitude is K = 
klT 2 V /M = h cS TV /h and t = t/T. 

In the case of linear coupling V — P. the echo-operator 
in Eq.(j8]) is therefore associated with the dimensionlcss 
unperturbed Hamiltonian 



H m = H KK + s m gp , 
and the perturbed effective Hamiltonian 



(30) 



(31) 



where e nm = g(s n - s m ) = gk T(s m - s„) is dimension- 
less. 

Replacing the quantum operators (0,p) with the clas- 
sical coordinates (6,p) in Hkr, one obtains the clas- 
sical Hamiltonian whose dynamics are governed by K. 
When K increases from zero the phase space struc- 
ture follows the Kolmogorov-Arnold-Moser (KAM) the- 
ory [39[ where the last invariant KAM torus is broken for 
K = Kr m 0, 97, and the motion becomes unbound for 
K > Kr. If K ~ 1, the classical phase space is mixed, 
and for K 3> 5 the classical motion may be considered 
completely chaotic having negligibly small stability is- 
lands 4l| . The classical counterpart of the Hamiltonian 
H, 



in Eq. ([30|) also presents essentially chaotic dynamics 
for K 3> 5 even when the strength of the linear cou- 
pling s m g is large, because the linear coupling simply 
represents a linear shift in kicked system's classical phase 
space. 

We verified the results in Sees. IIIII and llVl in two dif- 
ferent situations. The first is when the phase space of 
the kicked rotor is closed on the torus — tt < 6 < it, 
— 7r < p < tt so that the Hilbcrt space of the environment 
is finite (dimension N). Because the angle variable 8 
is bounded, the quantum momentum eigenstates are dis- 
crete, p\j) = h ctfJ \j) (j = -N/2, . . . , 0, A/2-1), and N is 
related to the effective Planck constant by H e g = 2tt/N. 
When the underlying classical dynamics of the unper- 
turbed Hamiltonians H m in Eq. (|3"0")) arc chaotic, at some 
finite relaxation time the environment's evolved state oc- 
cupies all the available momentum eigenstates. It is thus 
clear that in this case the effective size of the environ- 
ment's Hilbcrt space is N e g = N = 2ir/h c g, for all m. 
Substituting this into Eqs. (fT3| and (fT4|) one obtains 



1/A eff (x ftoff 



(Fr, 



(32) 



n=0 



The second situation is the case in which the envi- 
ronment Hilbert space is infinite (50l |. In this case the 
appearance of an effective dimension of the environment 
Hilbcrt space is well exemplified by the renowned phe- 
nomenon of dynamical localization[!(| EH]. It is well- 
known that during the time interval < t < t^, the 
kicked rotor Hamiltonian Hk with K > 5 and h c g « 
1 presents diffusion in the discrete momentum levels 
in good correspondence with the classical model, i.e., 
Ap 2 (t) fa D(K)t [where D(K) is the classical diffusion 
coefficient] |5lfl . But after a finite relaxation time scale 
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FIG. 1: (Color online) Typical behavior of the decoherence 
function F nm (t) as a function of time when the environment 
has chaotic classical dynamics. In these plots the environment 
corresponds to the kicked rotor on a cylinder with K = 5 
and linear coupling to the central system. The unperturbed 
Hamiltonian in the corresponding echo operator is given in 
Eq. (|30[l with s m g = 0.1 and the perturbed Hamiltonian is 
given by Eq. (|3ip with the strong perturbation e nm — 0.1. 
The time is measured in number of kicks. The full red line 
corresponds to h c fj = 0.92 and the dashed blue line corre- 
sponds to hes = 0.015. The initial state of the environment 
is w(0) = \p = 0)(p = 0|. 

tji, an important decrease in the diffusion rate is observed 
until the state ceases to spread in momentum space. This 
happens approximately when the occupation number in 
momentum space of the evolved state reaches the value 
iVeff = max„ m (iV m , N nm ) determined by the initial state 
w(0) (see Sec. Hi!]) , i.e., 

N cS ~ J w ^D(K)t R /h cff , (33) 

where J, called localization length, is essentially the 
width of the momentum eigenstates distribution. Due to 
the fully chaotic underlying classical dynamics, this phe- 
nomenon does not depend on the shape of the initial wave 
function (provided the initial state is not an eigenstate 
of the evolution operator). When dynamical localization 
takes place, the mean level spacing between quasi-energy 
eigenstates involved in the expansion of the localized (in 
momentum) initial wave function is A w lix j J . Thus, ac- 
cording to the Heisenberg principle, the minimum time 
required for the dynamics to resolve this level spacing is 
Ir~ 1/A oc J. Substituting tpt into Eq. (|3"3"]) one finally 
obtains 

l/N cS ~ 1/J cx h 2 cS => (F nm ), AF nm ex h 2 cS . (34) 

where we use the Eqs. ([T3|) and (fT4"]) . 

The behavior in time of the decoherence function 
Fnm(t), determined by the echo dynamics with the con- 
ditional unperturbed Hamiltonian given in Eq. (|30[) with 
K = 5, is illustrated in Figure [TJ in which we consider 
a strong coupling between the central system and the 
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FIG. 2: (Color online) Long-time behavior of the average 
{F n m(t)) (red solid lines) and the standard deviation AF nm 
(blue dashed lines), with n — m + 1, as a function of the 
coupling strength to the environment. The environment was 
modeled by a kicked rotator with K = 5, linearly coupled to 
a central system with s m — 1. The upper plot corresponds 
to the case where at large times the evolved state of the en- 
vironment spreads over all the available Hilbert space (the 
kicked rotor is defined on a torus) and in the bottom plot we 
simulate the kicked rotor on a cylinder, where dynamical lo- 
calization is observed. Other parameters are h e s = 0.025 (top 
plot) and fi c ff = 0.083 (bottom plot) and the initial state of 
the environment is cj(0) = \p = 0)(p = 0|. The three regions 
identified in the plots are explained in the text. 



environment. In both the case in which at large times 
the evolved state of the environment spreads over all the 
available Hilbert space and the case in which the available 
Hilbert space is much greater than the dynamical local- 
ization length, we found that at long-time scales, the de- 
coherence function F nm (t) fluctuates around an asymp- 
totic mean value (F nm (t)). This is typical in the long- 
time behavior of the fidelity amplitude in systems with 
fully classical chaotic dynamics in finite Hilbert spaces 
[27| . We confirmed that this behavior is independent of 
the initial state of the environment, w(0). 

The values of (F nm (t)) and AF nm for long times are 
plotted in FigJ2]as a function of coupling strength g. The 
top graph corresponds to the situation in which the en- 
vironment's phase space is closed on a torus and the 
bottom graph simulates the situation in which the en- 
vironment has an infinite Hilbert Space and dynamical 
localization is observed. The red solid lines correspond 
to the numerical results for (F nm (t)} and the blue dashed 
lines correspond to AF nm . From these plots it is possible 
to identify three echo perturbation regimes: region I is 
the weak perturbation regime and there is essentially no 
echo decay; Region II is an intermediate regime in which 
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FIG. 3: (Color online) (F nm ) (bottom graph) and AF nm (top 
graph) as a lunction of the effective Planck constant h e g for 
n — m + 1. The environment is a kicked rotator with K = 5 
with a strong linear coupling to the central system, e nm = 0.1. 
Blue circles correspond to a finite Hilbert space and the red 
filled dots correspond to the case in which dynamical localiza- 
tion is observed. The linear fits are drawn. The corresponding 
slopes of these fits are in good agreement with analytical pre- 
dictions and are given by: 2.03 ± 0.03 (top graph, filled red 
dots); 1.1 ±0.02 (top graph, blue circles); 2.03±0.03 (bottom 
graph, filled red dots); 1.08 ± 0.02 (top graph, blue circles). 



the mean value and standard deviation of the dccoher- 
ence function depend on both heti and the perturbation 
t m „ = g(s n — s m ); and region III is the strong pertur- 
bation regime in which the mean value and standard de- 
viation do not depend on the perturbation strength e nm . 
It is in this third region that the results in Eas.(|13p and 
(|14[) apply (with C = G for the particular initial state 
w(0) considered). In order to confirm these relations 
we plot in Fig|3]thc numerical calculation of log w (F nm ) 
(bottom graph) and log 10 (AF nm ) (top graph) as a func- 
tion of log 10 heff for the kicked rotor closed on a torus 
(blue circles) and in the case in which the Hilbert space 
is much larger than the localization length J (red filled 
dots). The linear fittings through the numerical points 
indeed confirm the relations given in Eq. (|32|) for the first 
case and the relations given in Eq.(|34p for the second 
case. Therefore, in the long-time regime, when the envi- 
ronment has fully chaotic underlying classical dynamics 
and the system-environment coupling is strong, all off- 
diagonal matrix elements of the reduced density matrix 
in Eq.([4"|) tend to zero in the semi-classical limit, i.e., 



CO 




FIG. 4: (Color online) Typical time evolution of the Hilbert- 
Schmidt distance Dhs between the evolved reduced state of 
a central system p(t) and the totally decohered state (p) in 
Eq. (|15[) for h c a — 0.53 (red upper curve) and fi e ff = 0.01 
(blue lower curve). In these graphs the central system is a 
harmonic oscillator and the system plus environment (kicked 
rotor) initial state is p(0) ®oj(0) where p(0) is a Schroedinger- 
cat-like state and uj(0) = \p = 0}(p = 0|. 
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FIG. 5: (Color online) Time-average value of the Hilbert- 
Schmidt distance Dhs between the evolved reduced state of 
a central system p(t) and the totally decohered state {p(t)) 
in Eq. (|15p . as a function of the effective Planck constant h e g . 
The central system is a harmonic oscillator and the environ- 
ment is a kicked rotator with K = 5 on a torus (blue circles) 
or on a "cylinder" (red filled dots), in which case dynamical 
localization is observed. The lines correspond to linear fittings 
of the numerical results whose slopes are given by 0.9 ±0.1 
(fitting over filled red dots) and 0.53 ± 0.01 (fitting over blue 
circles). The interaction between central system and envi- 
ronment is linear with e nm = 0.1 and the initial state of the 
system plus environment is the same as in Fig. [4] 



(F nm ), AF nm h ^°0. (35) 

In order to see the equilibration of the evolved reduced 
state of the central system p(t) to the totally decohered 
state /5 eq = (p(t)) in Eq. (fT5]) . we consider as a central sys- 
tem a quantum harmonic oscillator and an initial state 



p(0) = where \ip) = Af (\a) + \—a}) (\a) a coher- 

ent state, i.e, an eigenstate of the annihilation operator 
[29| ) . We consider this Schbdinger-cat-like state initially 
uncorrelated from the initial state of the environment 
w(0) = \p = 0){p = 0|. Equivalent results were obtained 
for different initial states p(0) <g> <£>(0). 
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The general behavior of the upper bound of the trace 
distance in Eq. (fTT]) is shown in Fig __] where we see that 
the Hilbert-Schmidt distance Dhs in long-time scales 
fluctuates around a mean value that decreases with fi, c g. 
In FigJS] we plot this mean value, given by the time- 
average (DHs{p(t),p B q)), as a function of the effective 
Planck constant. According to the result in Eqs. (f5^|) 
and (|34[) , the upper bound of the trace distance obtained 

_____ 1/2 

in Eq.(|19p has a dependence on hj s when the evolved 
state of the environment is spread over the entire avail- 
able Hilbert Space and has a dependence on h e g when 
the evolved state dynamically localizes before filling the 
available Hilbert space. This is confirmed by the linear 
fit over the numerical points in Fig. ([5]). Hence, in the 
semiclassical limit we obtain for the trace distance: 

(D(p(t),p cq )) h ^°0. (36) 

This means that in the semi-classical regime the central 
system always equilibrates to the totally decohered state 
given in Eq. fp~5|) , irrespective of the initial states of the 
system and the environment. 



VII. CONCLUSION 

We studied the decoherence process of a generic quan- 
tum central system coupled via a dcphasing-type inter- 
action to an environment whose underlying classical dy- 
namics is chaotic. We have shown that if the environment 
has chaotic classical dynamics and the coupling to the 
central system is strong, the time-average of the decoher- 
ence functions and the width of the fluctuations around 
this time-average are inversely proportional to the effec- 
tive Hilbert space dimension of the environment, defined 
as the dimension of the largest subspace onto which the 
projection of the initial state of the environment is not 
negligible, i.e., the maximum number of eigenstates of 
the perturbed or unperturbed conditional environment 
evolution operators needed to write the initial state of 
the environment. This means that if the effective Hilbert 
space of the environment is large, decoherence will oc- 
cur, i.e., in the long-time regime the coherences of the re- 
duced density matrix of the central system written in the 
preferential basis will decrease significantly and on aver- 
age will suffer only small fluctuations. If the quantum 
environment has chaotic underlying classical dynamics, 
the inverse of the its effective Hilbert space dimension 
will in general be proportional to some positive power 
of the effective Planck constant h e g, which implies that 
decoherence of the central system is guaranteed in the 
semiclassical regime (h ff << 1). This is true even for 
an environment with few degrees of freedom. We stress 
that, in general, if the coupling is not strong the central 
systems may have revivals of the coherences. 

We thus confirm the intuitive notion that the funda- 
mental quantity regarding the production of decoherence 
by chaotic environments is the dimension of the Hilbert 



space over which the environmental state is spanned, in- 
dependent of the number of environmental degrees of 
freedom Q. In [T__] the authors also discuss the con- 
nection between the effective Hilbert-space size and de- 
coherence, but in our case it becomes clear that the deco- 
herence is only guaranteed if the environment has chaotic 
underlying classical dynamics. 

For the system-environment models considered, in gen- 
eral the evolution of the reduced state of the system is 
not Markovian once the quantum information contained 
in the coherences of the reduced density matrix in the 
preferential basis flows back and forth between system 
and environment. However, we showed that the deco- 
herence process that occurs leads the central system to 
equilibrate in the semiclassical limit to a diagonal state 
given by the time-average of its evolved reduced density 
matrix. In this case, equilibration takes place indepen- 
dent of the initial state of the system and for a generic 
initial state of the environment (considered initially sep- 
arable from the s yste m) . This result is in agreement with 
that obtained in [___|. However, in our case, it is only the 
underlying chaotic classical dynamics of the environment 
that guarantees equilibration. 

We also investigate the entanglement decay of two 
nonintcracting central systems coupled via dcphasing- 
type interactions to identical quantum local environ- 
ments with classical chaotic dynamics. In this case we 
show that the reduced evolved state of the system equi- 
librates to a non-diagonal state in the preferential basis, 
which in the semiclassical limit (7i e __ — > 0) is arbitrarily 
close to a diagonal state (in the preferential basis). Thus, 
if the equilibrium state is not already disentangled, i.e., 
the evolved reduced state presents entanglement sudden 
death, a very small perturbation is sufficient to lead the 
system into a separable state. 

We confirmed all our analytical results with numerical 
simulations in which the environment is modeled by a 
kicked rotor in the regime where its underlying classical 
dynamics can be considered completely chaotic. 
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Appendix A: Derivation of the results in Eqs.(|13p 
and (fT4|) 

In this appendix we calculate the time-average and 
standard deviation of the decoherence functions F nm (t) 
assuming chaotic underlying classical dynamics of the en- 
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vironmcntal degrees of freedom. From here forward we 
assume n ^ m. We expand the conditional unperturbed 
and perturbed evolution operators, U m and U n in their 
own eigenbasis: 



fnm(t) = Tr [U*fi n w(p) 



(Al) 



ul = Y. e<lt/h \^)^ ( A2 ) 

Un = Y, e ~*' t/h $i)$l\> (^3) 
1 = 1 

where and are quasi-energies. To simplify notation, 
we do explicit the indexes n and m in 1-0;) and 

|^;). The sums in Eqs. (|A2p and (|A3[) may go to infinity 
if the spectrums of evolution operators are unbounded. 
We also expand the initial state of the environment in 
these eigenbases, as in Eqs. ([TU|) and (|lll) . Inserting (|A2[) 
and (|X3| in (O we obtain: 



1=1 l' = l 



where 



o w = {M$v) 



N„ 



(A4) 



(A5) 
(A6) 



k=l 



k' = l 



Observing Eq. (|A6[) it is possible to sec that matrix B^ 
is essentially non-zero in the intervals 1 < I < N m and 
1 < I' < N nm . We remind the reader that 1 < n,m < 
Nc, where Nc is the number of terms that significantly 
contribute to the expansion of the initial state of the 
central system p(0) in the eigenbasis {|«)}. Hence, as 
in Sec. we use the definition of effective dimension 
of the environment subspacc: N e s = max nrn (N m , N nm ) , 
where the maximum is taken over all possible n and m. 
Indeed, at any time, the reduced density matrix of the 
environment can be written as 



m = 



EE 

i=i k=i 



u lk (0)e-^-^ h \4>i)m = (A7) 



^^^ fc (0) e - i ^^^'/' i |^)(V fc | • (A8) 



1=1 k=l 

and the allegiance amplitude can be written as 



fnm(t) — EE 



-i(ii>-ii)t/n 



OivBv 



(A9) 



i=i i'=i 



For the time average of the decoherence functions we 
immediately obtain 



(A10) 



We consider a non-degenerate chaotic quasi-energy spec- 
tra and strong enough perturbation such that the bases 
{I'i/'i)} and {IV'i)} are uncorrelated. In this case, one can 
assume the overlap matrix O, formed by the elements 
Ow , to be a random unitary matrix [45| . In the limit of 
large N e g (such that (N c g + l)N c g sa A^ ff ), the matrix el- 
ements Ow can be considered complex random numbers 
with a Gaussian distribution oc exp(— N e g\Ou> | 2 ) [43| . 
We can then average over this distribution of matrix el- 
ements: 



\O u 



Oir_ = O*, = 0; 



\O w 



2/Nf 



off' 



(All) 



where • denotes the average over the Gaussian distri- 
bution of matrix elements. Using expressions (|A1 1|) in 
(|A10p . for TVcff > 1 we finally obtain: 



(-^77,777,) 



C 
iVef? 



(A12) 



where C = 7£ inv (w(0)) + Tr(w(0) 2 ) and K inv (ui(0)) = 
S£ZJ w i'(0)| 2 i s the generalized inverse participation ra- 
tio [271 ] of the initial reduced density matrix of the envi- 
ronment in the eigenbasis of the unperturbed evolution 
operator U m . 

In order to calculate AF nm in Eq.fJTJ]), we must calcu- 
late (F% m ) starting from: 



F. 



f,j,k,l=l f',j',k',l'-- 



.iitf-ti+tk-ti-Sf+tj'-ik'+Si')* 



x B* ff ,B j:jl Bl kl Bu,OffiO* 



'O k k'O w . 



(A13) 



For chaotic quasi-energy spectra it is reasonable to as- 
sume that they have non-degenerate gaps [? , 
so there exist only nine possible sets of conditions on 
{f 7 j,k,l} and {f',j',k',l'} for which the time average 
( e i(f/-ei+efc-Si-e f '+i'j'-€) e '+€!')*/ft) j s non-zero. In each 
case, one must compute the average over the Gaussian 
distribution of the elements of the overlap matrix O and 
then sum the results. In the limit of large N e g , we obtain, 



iV el 



>>+7v^E E 

JV cff k. 



Mo)| 2 | w (o) fefc | 2 + 



and 



therefore 



1 Z(^fc) = l 

Eq. (fT4|) where 



^ff 



(A14) 



v/2(^ inv (w(0))) 2 

E,E fc (Mo)| 2 - 



2Ez Mo)l 4 
K,(o)| 2 ). 



G 

and ip 



1=1 v = i 



Appendix B: The equilibrium state (p(t)} 

In this appendix we calculate the time average of the 
reduced system density matrix given in Eq.Q. The re- 
sult is the state in Eq. (fT5j) that was shown that corre- 
sponds to the equilibrium state of the central system in 



the semiclassical regime. We begin by expanding the con- 
ditional unperturbed operators and U n as in (|A2p and 
inserting these into Eq.(j4]) to obtain 

#*) = EE A ™ e-H( A ^)t(^ (0 )| Z )| m )( n | , (Bl) 

where A nm ; = e„ — sr m + Ej — E\ m ' and are the 

quasi-energies of the unperturbed conditional evolution 
operator U m . Taking the time average defined in Eq.© 
the only nonzero terms are those where A nm i = 0, that 
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